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Abstract: 

We propose a model to study the importance of double parton scattering (DPS) in 
Mueller-Navelet jets production at the LHC which is consistent with the BFKL frame¬ 
work used to compute the single parton scattering contribution to this process. We study 
this model in kinematics corresponding to existing and possible future measurements at 
the LHC and estimate the importance of this DPS contribution on relevant observables 
for this process, namely the cross section and the azimuthal correlation of the jets. 



1 Introduction 


The study of the production of two jets separated by a large interval of rapidity at 
hadron colliders was suggested by Mueller and Navelet [T] as a possible test of QCD in 
the high energy limit, which should be described by the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) approach |2, 3j 0J [5]. This pioneering study showed that, when the longitudinal 
momentum fractions of both jets are kept fixed, a BFKL calculation at leading logarithmic 
(LL) accuracy predicts a strong rise of the cross section when the rapidity separation 
between the jets increases, which is not expected in a fixed order approach. Later, another 
observable was proposed: the azimuthal correlation of the jets !M- The physical picture 
is that in a pure leading order collinear treatment the two jets would be emitted exactly 
back to back, while the multiple emissions of gluons in the rapidity interval between the 
jets in the BFKL framework should lead to a loss of angular memory and so a stronger 
decorrelation. 

Recently the azimuthal correlations of Mucller-Navelet jets were measured for the first 
time at the LHC by the CMS collaboration j8j|. It was found that this measurement can 
be described with reasonable accuracy [9] by a next-to-leading logarithmic (NLL) BFKL 
calculation supplemented by the use of the Brodsky-Lepage-Mackenzie procedure in 
adapted to BFKL dynamics du m, na E] to fix the renormalization scale (see also m 
m mi hb]). However, at high energies and low transverse momenta, which is the region 
where BFKL effects are expected to be enhanced, parton densities can become large 
enough that contributions where several partons from the same incoming hadron take 
part in the interaction could become important. These multipartonic interactions (MPI) 
would make the comparison of BFKL calculations with experimental data much more 
difficult. The goal of the present paper is to propose a model to study MPI effects 
in this process and to evaluate their importance in the kinematics of the CMS study 
and of potential future studies at the LHC. We will restrict ourselves to the case where 
there are at most two scattering subprocesses and where both these scatterings are hard. 
This is generally known in the literature as double parton scattering (DPS). We note 
that some work was already done in this direction, where it was found that the DPS 
contribution could become more important than the single parton scattering (SPS) one 
in some kinematics [T9]. However the DPS contribution was evaluated in the collinear 
framework which makes the comparison with a BFKL calculation questionable. Here we 
propose a model to evaluate the DPS contribution in a way which is consistent with the 
BFKL framework. 


2 Formalism 

2.1 BFKL calculation 

A detailed description of the BFKL calculation used to compute the SPS contribution 
to Mueller-Navelet jets production can be found in Refs. [29, HU 9J H2J- For similar 
treatment see [23, 25 HU 26] 
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Here we will only illustrate a few key points that will serve as a basis to build our 
DPS model. 

A diagrammatic representation of this process is shown on Fig. [lj Both incoming 
protons radiate a collinear parton which couples to the BFKL Green’s function G, de¬ 
scribing the multiple emission of untagged gluons, and the corresponding outgoing jet via 
the jet vertex. The differential cross section for jets with transverse momenta kjy and 
rapidities yj i% reads 


da 

d | kyi | d | k j >2 1 dy yi dyj t 2 


h u 


d 2 k 1 d 2 k 2 $(x.7 ) i, -ki)G(ki, k 2 , s)<f>(a:j !2 , k 2 ), 


(1) 

where the impact factors $ are a convolution of the jet vertex V with the parton distri¬ 
bution function (PDF) /: 


kj) 


J dxif(xi)V(ki,Xi). 


( 2 ) 


The leading order jet vertex is 


H a { 0 ) ( k i,Xi) = (l - ) |kj,i|<5 (2) (ki - kjy), 


where Ca = iV c = 3 is to be used when the proton radiates a gluon and Cp 
1)/(2 N c ) = 4/3 when it radiates a quark. 


(3) 





Figure 1: Kinematics of the SPS contribution in the BFKL approach. 


2.2 DPS model 

To model the DPS contribution to this process, we want to stay as close as possible to 
the spirit of the BFKL approach described in the previous section to avoid discrepancies 
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which could arise due to using different frameworks. We illustrate the DPS contribution 
on Fig. [2j In this model, each jet is emitted by a different BFKL-like ladder. Each 
incoming proton emits two partons: one collinear, which enters the jet vertex “close” to 
its fragmentation region, and a gluon with nonzero transverse momentum which initiates 
the ladder attached to the other proton. 



In practice, to compute this DPS contribution, one would need to introduce some 
kind of “hybrid” double parton distributions, related to the probability to emit both a 
collinear parton and a gluon with some transverse momentum. Since at the moment 
almost nothing is known about such distributions, we use instead the simple factorized 
ansatz to compute the DPS contribution according to 

^DPS = - j (4) 

0eff 

where <Xf w d(bwd) is the inclusive cross section for one jet in the forward (backward) direc¬ 
tion and cr e fr is a phenomenological quantity related to the density of the proton in the 
transverse plane. According to measurements at the Tevatron [27U28U29U30] and at the 
LHC [3T| [32], cr eff should be of the order of 15 mb but there is some discrepancy between 
these measurements. To account for this uncertainty we will vary <j e g between 10 and 
20 mb in our calculation (note that this already amounts to an uncertainty of a factor 
2 on the DPS cross section). By using the approximation in Eq. Q, we have reduced 
the problem to the calculation of the cross section for the production of one jet in the 
forward (or backward) direction. 

To compute this cross section we start from the BFKL cross section for dijet produc¬ 
tion in Eq. ([!]) in the LL approximation, which is illustrated on Fig. [3] (L). Our goal is 
to compute the cross section for one forward jet, as shown on Fig. [3] (R). For this we 
replace the part of the diagram containing the coupling to the lower proton, the lower 
jet vertex and the Green’s function by an unintegrated gluon distribution (UGD), which 
couples directly to the lower proton and describes the multiple gluon emission between 
this proton and the jet. For simplicity we use the LO jet vertex as written in Eq. (J3|, 
therefore the transverse momentum probed in the UGD is the one of the produced jet 
and its longitudinal momentum fraction value x is given by the condition for the outgoing 
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Figure 3: Left: Muellcr-Navelet jets production at LL accuracy. Right: Inclusive forward 
jet production. 

parton, which will then form the jet, to be on-shell. Unintegrated gluon distributions are 
much less known than usual collinear PDFs. Several models have been proposed, see for 
example Refs. |33, [33j [35], [36], 3T] [38, [39] HD] , but they can differ significantly in some 
phase space regions. In particular, we observe that for inclusive forward jet production 
the normalization depends strongly on the parametrization chosen. For this reason we 
write the differential cross section in kx —factorization [01], 02, T3[. 00, 05| 06, 37| as 



(5) 


where the unintegrated gluon distribution is normalized according to 



( 6 ) 


and K is a dimensionless normalization factor which we will fix independently for each 
UGD parametrization with LHC data as explained in the following. The cross section 
for one forward jet production as a function of the transverse momentum of the jet was 
measured by the CMS collaboration at a center of mass energy of 7 TeV in the range 
3.2 < \yj\ < 4.7 [08j. For each parametrization, we determine the range of K compatible 
with the CMS measurement in the lowest transverse momentum bin. We choose to fit 
data only in the first bin because none of the parametrizations can describe perfectly the 
data over the whole px range and in the following we will be interested in rather low 
transverse momenta. The comparison with CMS data after fixing K according to this 
procedure is shown on Fig. [4] for several choices of UGD parametrizations: KMR [35] . 
AO [38], KS [39] and JH2013 [40]. We see that all models describe reasonably well the 
trend of the data. 

3 Results 

Now that we have all the necessary ingredients to compute both the SPS and the DPS 
contributions to this process, we are ready to estimate their relative importance (in the 
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Figure 4: Comparison of CMS data at y/s = 7 TeV for inclusive forward jet cross section 
in the range 3.2 < yj <4.7 with estimates based on several UGD parametrizations (see 
text for references) after fixing the parameter K. 


following we will neglect any possible interference between them). First we will compare 
the magnitude of the cross sections for both contributions. Then we will study how the 
DPS contribution could affect the azimuthal correlations of the jets. We will focus on 
four choices of kinematical cuts: 

• y/s = 7 TeV, k/^ = |kj )2 | = 35 GeV, 

• \/i = 14 TeV, |k,/i| = |k/ j2 | = 35 GeV, 

• \/s = 14 TeV, |k/j| = |kj )2 | = 20 GeV, 

• Vs = 14 TeV, |k Ja | = |k J)2 | = 10 GeV. 

The first choice is similar to the cuts used by the CMS analysis of azimuthal correlations 

of Mucller-Navelet jets at the LHC jS]. It is important to evaluate the DPS contribution 
for these kinematics to make sure that the good agreement found in Ref. pj is not just 
a coincidence. The other three choices correspond to the higher center of mass energy 
that the LHC is expected to reach soon. The choice |kyi| = |kj )2 | = 35 GeV allows us 
to compare with the results at 7 TeV to evaluate the importance of a change of y/s. The 
last two choices correspond to lower transverse momenta at which measurements could 
become possible in the future. They are particularly relevant since MPI are expected to 
become more and more important at lower transverse momenta. The rapidities of the jets 
are restricted according to 0 < yy\ < 4.7 and —4.7 < yj ;2 < 0. We use the MSTW 2008 
parametrization |49j for collinear parton densities. To estimate the uncertainty associated 
with the choice of the UGD parametrization needed to compute the DPS cross section, 
we use the same four parametrizations as in Fig. [4] and take the extreme values as an error 
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Figure 5: Comparison of the differential cross section obtained at LL (green) and NLL 
(red) accuracy in the BFKL approach and the DPS cross section (blue) for the four 
kinematical cuts described in the text. 


band. In the case of the NLL BFKL calculation, for which a jet can be made of more 
than one parton, we use the anti -k t jet algorithm [50] with a size parameter R = 0.5. 

3.1 Cross section 

In this section we compare the cross section obtained with the DPS model described in the 
previous section with the SPS contribution computed in the BFKL approach. On Fig. [5] 
we show, for the four kinematical cuts described above, the SPS cross section computed 
in the BFKL approach, either in a pure LL treatment or when including NLL corrections 
both to the Green’s function [SB E2] and to the jet vertex [SB EU ESI ESI EZ] ; together 
with the DPS cross section. The first observation is that the uncertainty is much larger 
for the DPS cross section than for the SPS one. This uncertainty is due to the variation of 
cr e ff between 10 and 20 mb and the spread between the different UGD parametrizations. 
The cross section obtained in the three treatments increases at larger energy and smaller 
transverse momenta, especially for large rapidity separations. This is important because 
this means that a new analysis in the spirit of Ref. [8j could provide much more statistics 
with the same luminosity. The DPS cross section is always smaller than both the LL 
and NLL BFKL calculations, but this difference reduces with increasing Y and smaller 
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Figure 6: Ratio of the DPS cross section to the LL (green) and NLL (red) BFKL SPS 
cross section as a function of Y for the four kinematical cuts described in the text. 

transverse momenta. This can be seen more clearly on Fig. [6] where we show the ratio of 
the DPS cross section to the SPS one. This ratio grows faster with Y when one treats the 
SPS contribution at NLL accuracy than at LL accuracy. This is due to the well-known 
fact that LL BFKL predicts a very strong rise of the partonic cross section with increasing 
rapidity separation between the jets. Note however that this LL calculation is not reliable 
anyway since for example it cannot describe CMS data on azimuthal correlations of the 
jets i but we show it for illustrative purposes. It is important to note that for the first 
choice of kinematical cuts the DPS contribution is always smaller than the SPS one by 
at least one order of magnitude. As a consequence, the comparison made in Ref. [9] 
with CMS data should be safe. The situation could be different for other kinematical 
conditions. For example, for y/s = 14 TeV, | k^i | = |kj i2 | = 10 GeV, at the edge of the 
uncertainty band the DPS contribution can become almost as large as the NLL BFKL 
cross section. However the uncertainty is again very large so this ratio is also compatible 
with much smaller values of the order of 1%. 

3.2 Azimuthal correlations 

Even if the DPS cross section is always smaller than the SPS one in the kinematics under 
study here, there could be a measurable effect of DPS on more exclusive observables like 
the azimuthal correlation between the jets. Since our simple DPS model neglects any 
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correlation between the jets, the inclusion of this contribution can only lead to a larger 
decorrelation than what is found in an SPS study, the size of this effect depending on the 
magnitude of the DPS cross section compared to the SPS one. 

To measure the azimuthal correlation between the jets, one can define the angular 
coefficients (cos rup), with ip = (j)jj — (ftj j2 — 7r. In particular, n — 1 corresponds to 
(costp) for which a value of 1 corresponds to always back-to-back jets while a value of 0 
corresponds to completely uncorrelated jets (as we have for the DPS contribution). The 
correlation for all (SPS and DPS) events is given by 


(cos rup) = 


osps (cos ntp) sps + 0 - DP s(cos?r<p) D ps 


T3PS + cr D ps 


(cos rup) 


SPS 


| °DPS 
O'SPS 


(7) 


where it is clear that 

0 < (cos rup) < (cosn<p)sps • (8) 

On Fig. [7] we show the first of these correlation coefficients, (cos<p), with and without the 
inclusion of DPS effects, for the four choices of kinematics described previously. I 11 all 
cases the SPS and SPS+DPS calculations are compatible with each other. The reason is 
that, as shown on Fig. [6j for some choices of the parameters the DPS cross section can 
be negligible compared to the SPS one. Even at the lower edge of the uncertainty band 
the deviation to SPS is not significantly larger than the uncertainty on the NLL BFKL 
calculation, except for the last choice of kinematics (y/s = 14 TeV, |kyi| = |k^ 2 1 = 10 
GeV) at large rapidity separations. 

One can also study higher harmonics (cos rup) with n > 1. This allows us to compute 
ratios of the kind (cos 2<p )/(cos <p) which have been found to be much less dependent 
on the choice of the scales than individual coefficients (cos rup), both at LL and NLL 
accuracy [5HI EH EDI EJ EH El], and to show a good agreement with experimental data 
at NLL accuracy, even without using any scale-fixing procedure [9]. These ratios are not 
affected by the uncorrelated DPS contribution under study here since, as can be seen 
from Eq. (jTj) , the quantity 1 + will disappear in any such ratio. 

Another observable related to the angular correlation of the jets is the azimuthal 
distribution yyy, which can be computed making use of the previously defined coefficients 
(cos rup) according to 


-Sr = -i 1 + 2 E 

cr do? 2n I ^ 

' l n= 1 


cos (rup) (cos (rup)) 


(9) 


Since in our model the DPS contribution produces uncorrelated jets, the corresponding 
azimuthal distribution is a constant. Similarly to the case of (cos ip) studied above, 
the total azimuthal distribution can be obtained by taking the average of the azimuthal 






(cos ip ) 


(cos ip ) 



(cos ip) 


(cos ip) 




Figure 7: (cos<p) as a function of Y without (red) and with (blue) the DPS contribution 
for the four kinematical cuts described in the text. 


distribution for the SPS and DPS contributions weighted by their respective cross section: 

dasps dcrpps 
d <p dtp 

d^PS + ^DPs\ (1Q) 

^sps + ^dps \ dtp 2vr ) 

In Fig. [8] we compare the azimuthal distribution obtained by taking into account both 
the SPS and DPS contributions with the results of the SPS calculation alone in the 
range 8 < Y < 9.4. We observe that for the first three choices of kinematical cuts the 
maximal possible deviation from SPS induced by the DPS contribution is smaller than 
the uncertainty on the SPS calculation itself. On the contrary, for the choice a/s = 14 
TeV, |kjj| = |kj i2 | = 10 GeV, the DPS contribution can lead to a significantly flatter 
distribution but there is an important overlap between the two treatments which are 
therefore still compatible with each other. 


1 da 1 

a dip osps + ctdps 


3.3 Determination of the DPS fraction from experimental data 

In the previous sections we discussed the importance of double parton scattering for sev¬ 
eral observables in this process. However we found that our estimation is affected by a 
large uncertainty, which does not allow us to reach firm conclusions. To estimate the 
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Figure 8: Azimuthal distribution integrated in the range 8 <Y< 9.4 without (red) 
and with (blue) the DPS contribution, for the four kinematical cuts described in the text. 


actual DPS contribution, one could try to extract the DPS fraction from future experi¬ 
mental data under the assumption that the DPS contribution is fully uncorrelated. As 
already explained above, such a DPS contribution would not affect conformal ratios like 
(cos2(p)/(cos<p). Therefore, if a future measurement shows that the BFKL calculation 
describes this observable well but overestimates significantly (cos ip ), this could be inter¬ 
preted as a large DPS fraction, which could be estimated by inverting the relation in 
Eq. 0: 

^DPS (cos(p)sps 


^SPS 


(cos <p) 


- 1 , 


( 11 ) 


where (cos<p)sps is the value obtained with the BFKL calculation and (cos <p) is the value 
determined experimentally. 
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4 Conclusions 


In this paper we proposed and studied a model to estimate the importance of double 
parton scattering in the production of two jets separated by a large interval of rapidity 
at hadron colliders. This simple model is based on the simple factorized ansatz which 
neglects any correlation between the jets emitted by two separate ladders. We esti¬ 
mated the uncertainty of our calculation by using different unintegrated gluon density 
parametrizations and by varying the value of a e s in a range compatible with experimental 
determinations of this quantity. It turns out that the resulting uncertainty on the DPS 
cross section is rather large. Still, this cross section is always smaller than the SPS one 
in the LHC kinematics we considered here. We also studied the impact of double parton 
scattering on the angular correlation between the jets. We found that the inclusion of this 
contribution leads to predictions still compatible with a NLL BFKL calculation including 
only single parton scattering. However, if one considers the set of parameters giving the 
largest DPS contribution, for low transverse momenta and large rapidity separations the 
effect of DPS can become larger than the uncertainty on the NLL BFKL calculation. 
Therefore in this region a more careful analysis or experimental data would be required 
to conclude. Finally, since our study relies on the relative size of the SPS and DPS 
cross sections, we would like to stress that a measurement of the absolute cross section 
for Mucllcr-Navelet jets production in kinematics similar to [8j, where DPS effects are 
expected to be negligible, would be extremely valuable to confirm these results since it 
would provide a check of the SPS cross section normalization. 
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